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ABSTRACT 

A  technique  is  presented  for  predicting  the  response  to  failure  of 
composite  laminates.  In  this  technique,  cubic  spline  interpolation 
functions  are  used  to  represent  basic  stress-strain  curves  obtained  from 
simple  tests.  This  representation  in  conjunction  with  an  incremental 
constitutive  law  is  used  to  generate  response  of  laminates  subjected  to 
incremental  loads.  The  ultimate  load-carrying  capacity  of  the  laminates 
is  determined  by  the  plywise  application  of  a  new  failure  criterion 
using  strain  energy  under  longitudinal,  transverse,  and  shear  loadings 
as  independent  parameters.  Results  obtained  by  this  technique  are 
compared  with  the  available  experimental  data  of  boron-epoxy  laminates 
as  well  as  the  responses  generated  by  methods  of  Petit-Waddoups  and  Hahn. 
Analytical  results  determined  by  the  present  theory  correlate  with  the 
experimental  data  better  than  those  obtained  by  the  other  two  techniques. 
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SECTION  I 
INTRODUCTION 

Filamentary  composite  laminates  are  being  used  in  aircraft  structural 
components  in  increasing  numbers.  To  use  these  materials  satisfactorily 
we  need  to  know  their  stress-strain  behavior  and  ultimate  strength  under 
design  loads.  This  data  can  be  obtained  experimentally,  but  this  is  both 
expensive  and  inconvenient,  since  a  large  variety  of  laminates  can  be 
fabricated  from  basic  unidirectional  plies.  The  other  technique  would 
be  to  analytically  relate  the  properties  of  the  lamina  to  those  of  the 
laminate.  This  would  require  knowing  the  material  biaxial  properties 
of  the  lamira  for  various  load  combinations,  which  information  is  not 
available  at  present. 

Presently  properties  of  the  lamina  under  simple  loading  conditions 
(uniaxial  tension,  uniaxial  compression,  and  pure  shear)  may  either  be 
estimated  from  the  properties  of  its  constituents  (Reference  1 )  or  deter¬ 
mined  experimentally.  The  estimation  procedure  has  not  proved  effective, 
so  mechanical  properties  of  the  lamina  are  obtained  by  tests.  The 
behavior  of  the  lamina  under  complex  stress  states  is  predicted  by  using 
some  criterion  to  relate  the  parameters  of  strength  obtained  from  simple 
tests  to  the  stress  states.  In  this  technique  we  normally  assume  that 
the  yield  stress  and  the  ultimate  strength  are  synonymous.  Some  of 
these  failure  theories  (Reference  2)  assume  an  interaction  between  the 
failure  modes  to  generate  a  smooth  failure  surface.  Others  treat  the 
failure  modes  as  being  independent. 

Even  if  the  behavior  of  the  lamina  can  be  accurately  predicted,  the 
problem  of  determining  the  performance  of  laminas  in  a  laminate  still 
remains.  The  degradation  of  strength  of  individual  laminas  may  not 
precipitate  failure  of  the  laminate.  For  example,  in  the  case  of  a 
cross  ply  (0790°)  laminate  subjected  to  normal  stresses  (o7,o2)  acting 
in  the  direction  of  the  fibers,  the  degradation  of  strength  of  laminas 
in  the  transverse  direction  does  not  result  in  the  failure  of  the 
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laminate;  it  is  capable  of  sustaining  stresses  higher  than  those  causing 
transverse  degradation  of  the  individual  laminas. 

This  lamina-laminate  response  has  been  a  subject  of  many  theoretical 
and  experimental  studies.  Tsai  (Reference  3)  assumed  the  lamina  to  be 
quasi -homogeneous,  orthotropic,  and  linearly  elastic  to  failure,  and  he 
determined  stress  states  in  all  plies  due  to  applied  loads.  To  satisfy 
Hill's  failure  criterion,  the  plies  were  considered  to  be  degraded  and 
their  stiffnesses  modified;  however,  in  this  technique,  the  exact  nature 
of  degradation  is  not  explicitly  defined.  Later  Tsai,  et  al  (Reference  4) 
to  account  for  post  yielding  behavior  of  the  laminate,  proposed  a  new 
technique  incorporating  some  of  the  assumptions  of  netting  analysis. 

After  initial  degradation,  this  technique  postulates  that  the  laminate 
can  sustain  stresses  in  the  direction  of  fibers,  which  is  made  possible 
by  the  assumption  that  an  internal  agency  exists.  This  agency  would 
provide  resistance  to  the  externally  applied  loads  by  developing  axial 
loads  along  the  fibers. 

Noyes  and  Jones  (Reference  5)  extended  the  concepts  of  Reference  3 
by  modifying  Hill's  criterion  to  determine  the  onset  of  failure  in  a  ply 
or  plies.  To  determine  the  associated  failure  mode,  they  proposed  an 
empirical  technique  based  upon  inequalities  involving  the  stress  state 
and  parameters  of  strength.  If  the  analysis  indicates  a  transverse  or 
a  shear  failure  mode,  both  the  transverse  and  the  shear  stiffnesses  of 
the  affected  plies  are  set  to  zero;  with  longitudinal  failure,  all  the 
stiffnesses  are  reduced  to  zero.  The  loading  is  continued  with  the 
remaining  plies  Intact  until  general  failure  of  the  laminate.  This  pro¬ 
cedure  assumes  instantaneous  failure  of  the  plies;  however,  the  failure 
may  be  gradual  with  loading  rather  than  Instantaneous.  Chiu  (Reference  6) 
suggested  the  approach  by  Noyes  and  Jones  be  modified  to  reduce  the 
affected  stiffnesses  and  strengths  gradually. 

All  the  above  approaches  assumed  linearity  of  lamina  stress-strain 
relations.  Petit  and  Waddoups  (Reference  7)  relaxed  this  restriction; 
to  determine  response  of  the  laminates,  they  used  properties  of  unidirec¬ 
tional  laminae  in  conjunction  with  a  piece-wise  linear  approach. 
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The  loading  is  proportional  and  is  applied  in  small  increments  to  compute 
corresponding  strain  increments.  These  increments  are  added  to  the 
previous  strains  to  obtain  the  resultant  strains  at  a  particular  load 
level.  These  strains  are  used  to  determine  the  elastic  constants  to  be 
used  in  the  next  cycle.  When  strain  components  reach  limiting  values, 
lamina  or  laminae  are  assumed  to  unload  rapidly,  thereby  transferring 
their  share  of  load  to  the  Intact  layers.  The  laminate  fails  when  it 
cannot  sustain  any  more  load  increments.  Tangent  moduli  are  obtained 
from  stress-strain  curves  using  the  lamina  strains:,  however,  this  does 
not  consider  the  fact  that  the  strains  in  the  input  stress-strain  curves 
are  simple,  while  the  lamina  strains  in  the  laminate  are  not.  (This 
aspect  will  be  discussed  in  greater  detail  In  Section  II,)  Additionally, 
stress-strain  curves  of  laminates  obtained  by  this  method  are  very  much 
dependent  upon  the  size  of  the  Increment  used  in  the  applied  loads. 

Recently  a  new  technique  (References  8  and  9)  has  been  proposed  for 
the  study  of  the  nonlinear  behavior  of  laminated  composites.  This  tech¬ 
nique  Incorporates  only  nonlinearity  of  shear,  represented  by  the  equation 

<6  s  S66  <r6  +  S6666  <r63  (!) 

where  ag,eg  are  shear  stress  and  shear  strain,  respectively,  and  Sgg,Sgg6g 

arc  constants.  This  approach  does  not  distinguish  between  tension  and 
compression  behaviors,  however,  and  It  does  not  deal  with  the  problem  of 
failure. 

This  report  describes  a  technique  using  cubic  spline  interpolation 
functions  to  represent  the  basic  stress-strain  data,  namely, 

a.  Longitudinal  tension  and  compression, 

b.  Transverse  tension  and  compression,  and 

c.  Shear. 

This  functional  form  of  stress-strain  curves  provides  an  accurate  repre¬ 
sentation  of  stresses,  strains,  and  moduli  over  the  entire  range  of  »he 
curves.  The  Incremental  constitutive  law  relates  the  Increments  of 
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stresses  and  strains.  The  basic  data  and  the  constitutive  relations  are 
used  to  determine  the  ultimate  cumulative  response  of  the  laminates  under 
incremental  loading  by  a  new  failure  criterion.  This  criterion  is  based 
upon  the  concept  that  strain  energies  under  longitudinal,  transverse,  and 
shear  loadings  are  independent  parameters.  The  analytical  response 
obtained  by  the  three  methods  (the  proposed  technique  and  those  of 
References  7  and  8  and  9)  is  compared  with  the  available  boron-epoxy 
experimental  results. 
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SECTION  II 
METHOD  OF  ANALYSIS 

The  basic  constituent  of  a  laminate  is  a  unidirectional  lamina  or 
ply.  This  lamina  is  essentially  orthotropic.  The  mechanical  properties 
of  the  lamina  under  static  loads  are  determined  by  conducting  five  simple 
tests  (uniaxial  tension  and  compression  along  and  transverse  to  the 
material  axes  and  shear).  Generally,  these  tests-indicate  a  nonlinear 
behavior  of  the  lamina.  For  this  reason,  an  incremental  type  constitutive 
law  is  desirable  to  define  the  response  of  the  lamina  under  general  states 
of  stress.  To  define  the  general  incremental  constitutive  relationship, 
it  is  assumed  that 

a.  The  increment  of  strain  depends  upon  the  strain  state  and  the 
increment  of  stress;  and 

b.  The  Increment  of  strain  is  proportional  to  the  increment  of  stress. 
Using  the  above  assumptions,  the  incremental  constitutive  law  can  be 
written  as 

d*l)  =  S«r»,*u)  d<Tr,  =  *.a.3»  (2) 

where  de^,  darJ  are  strain  and  stress  Increments  and  .  Is  a  function 
of  the  current  strains,  e^..  J 

The  Incremental  constitutive  relation  using  contracted  notation  for 
the  lamina  under  generalized  plane  stress  reduces  to 

dc,  *  S|J(«,)d<rJ  (lfJs  1,2,6)  (3) 
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Equation  3  written  explicitly  becomes 


d€, 

Sll  $12  S,6 

drr, 

d<2 

s 

S2|S22  s26 

d<r2 

d€(?1 

S6IS62  S66 

d<r6 

or  simply 

[d«]  =  [  s  ][d<r]  (5) 


where  dcr^ ,  de1 
do2>  de2 

do6’  de6 
S 1 1  *  ^i2* • 


normal  stress  and  strain  increments  in  the  fiber 
direction, 

normal  stress  and  strain  increments  in  the  transverse 
direction, 

shear  stress  and  strain  increments 

elements  of  compliance  matrix  representing  the 
average  values  during  the  increment  of  stresses. 


Since  the  lamina  is  assumed  to  remain  orthotropic  at  all  load  levels, 
the  elements 


u> 

Gi 

II 

^28  =  S6j  =  S62  s  0 

(6) 

SM 

= 

* 

»  ®22  =  *^E22  *  ^66  s  l^|2 

SI2 

s  s2l 

=  “  V\z/^\\  ~  “  *2l/E22 

(7) 

Combining  Equations  4,  6,  and  7  provides  incremental  stress-strain 
relations  for  the  lamina 


dc,  *  do-,11  -  v)2B)/e„ 


(8) 


6 


- . .  -■ ' 
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df2  =  d<r2(l 


where 


<*««*  d^«/6l2 


B  =  do*  /da; 


Equations  8  and  9  reveal  that  it  would  be  erroneous  to  use  de^  or 
de2  lepresenting  the  biaxial  stress  state  do-j  ,do2  to  determine  E^  or 
E22  from  stress  strain  curves  obtained  under  simple  loading  conditions. 
Fnr  example,  de2  of  Equation  9  corresponds  to  the  curve  ON  (Figure  1) 
on  the  plane  OEHG  while  the  simple  stress  strain  curve  OM  lies  on  the 
plane  OEDC.  Since  stress-strain  data  similar  to  ON  is  not  available, 
we  have  assumed  that  simple  equivalent  strain  increments  can  be  computed 
from  the  following  expressions: 


In  the  application  of  the  Incremental  constitutive  law  to  multi¬ 
directional  laminates,  Equation  5  for  the  k-th  ply  is  written  as 
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where 


Stiffness  matrix  of  k-th  ply 


Stress  increment  in  the  k-th  ply  relative  to  the 
material  axes  1,2 


Strain  increment  in  the  k-th  ply  relative  to 
the  material  axes  1,2 


In  the  case  of  multidirectional  laminates,  load  axes  X,  Y  generally  do 
not  coincide  with  the  material  axes  1,  2  (Figure  2).  Stress  and  strain 
increments  in  two  coordinate  systems  are  related  by  a  transformation 
matrix  Jt  jk>  i.e., 

os) 

C*1*]/  [T]k[d«]k  06) 

where  [dajk  and  [dejk  are  the  stress  and  strain  increments  in  the 
X,  Y  coordinate  system.  Substitution  of  Equations  15  and  16  in  14, 
yields 


HMkH 

Khhh 

[<  ]„H 


(17) 


. — — 
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il 


where 


Assuming  that  stresses  are  uniformly  distributed  through  the  thickness  of 
each  ply,  stress  resultant  increments  [dN J  ,  in  the  X,  Y  coordinate  system 
are  given  by  J 

M*  „?, ["*],'*  09) 

tk  -  thickness  of  k-th  ply 
P  =  number  of  plies  in  a  laminate 

On  substituting  Equation  17  in  Equation  19,  stress  resultant  increments 
become 


M-j'jnM, 


k1-  k 


Assuming  that  the  strain  Increments  [del,  are  the  sames  for  all  plies 
of  the  laminate,  i  .e. ,  '  *“ 

[d.]  =  [dl]  (2, 


Equation  20  becomes 


where 


[«]  =  [*]  [d.*] 
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Inversion  of  Equation  22  yields 

[dl]  =  [A][dN] 


(24) 


In  Equation  24  [a]  represents  the  average  compliance  properties 
of  the  laminate  during  the  (n+l)th  load  increment.  However,  [a]  is 
not  known  when  (n+1  )th  load  increment  is  applied.  To  overcome  this 
difficulty  elastic  properties  at  the  end  of  the  nth  load  increment  are 
used  in  Equation  24,  i.e., 

[^L,  ■[*!  M„+l  (25) 

The  strain  increments  [dl]n+,  obtained  from  Equation  25  are  used  in 
Equations  17,  15,  16,  12,  and  13  to  calculate  [dol,  [del  de  I 
and  deg |  These  stress  and  strain  increments  are  added  to  'Eq- 

stresses  and  strains  at  the  nth  load  increment  to  obtain  the  current 
stresses  and  strains  in  all  plies.  The  current  stresses  and  strains  are 
employed  to  determine  the  average  elastic  properties  of  the  plies  and  a 
new  [Aj  is  computed.  This  procedure  Is  repeated  till  the  difference 
between  two  values  of  [dljn+1  in  two  consecutive  cycles  is  small  to  the 
order  of  approximation  desired.  The  Equation  25,  J hen  becomes 

ML.  •[*!♦■  M/i,  '  <26> 

* 

The  repetitive  us  of  the  procedure  outlined  abr/ve  generates  the  stress- 
strain  response  of  multidirectional  laminates./ 
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SECTION  III 
FAILURE  CRITERTON 

The  incremental  loading  procedure  outlined  in  Section  I, to  predict 
the  response  of  the  laminate,  cannot  proceed  indefinitely.  A  stage  is 
reached  where  a  lamina  or  laminas  can  no  longer  sustain  additional  loads. 
The  failure  state  of  the  lamina  under  general  stress  states  is  determined 
by  a  criterion  relating  its  behavior  to  that  under  simple  load  conditions. 
In  the  formulation  of  the  criterion  for  anisotropic  materials,  essentially 
two  approaches  (Reference  2)  have  been  used.  One  approach  characterizes 
the  failure  state  to  be  when  stress/strain  components  exceed  the  limit 
obtained  by  simple  tests.  The  second  approach  has  been  to  use  different 
quadratic  forms  of  stresses  with  or  without  linear  terms  to  establish  the 
failure  state.  The  coefficients  of  the  quadratic  expression  are  evaluated 
by  using  results  of  simple  or  both  simple  and  biaxial  tests.  These 
approaches  assume  that  yield  and  ultimate  states  are  analogous.  For 
materials  exhibiting  nonlinear  response  under  loading,  the  two  states  are 
not  the  same,  and  the  failure  state  would  be  a  function  of  both  stress 
and  strain  states. 

A  scalar  function,  f,  defining  the  failure  condition  of  materials 
exhibiting  nonlinear  behavior  can  be  written  as 

f  (cr,  c,K  )  =1  (27) 

where  a,  e  are  the  stress  strain  states  and  K  is  the  material  character¬ 
istics. 

An  explicit  form  of  Equation  27  proposed  in  this  report  uses  the 
scalar  strain  energy  as  a  measure  to  determine  the  level  of  the  effect 
of  both  stress  and  strain  states  on  the  material  behavior.  Assuming 
strain  energies  as  independent  parameters,  the  failure  criterion  for 
orthotropic  materials  may  be  expressed  as 

Kl|  [I-  V*„]  . I  ■'.*.*>  (28) 

I) 
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where  £ij  are  the  current  strain  components  and  m  is  a  parameter  defining 
the  shape  of  the  failure  surface  in  the  strain  energy  space.  The  failure 
criterion  (Equation  28)  is  based  upon  total  strain  energies.  It  includes 
the  effect  of  hydrostatic  loading.  This  inclusion,  according  to  Reference 
10,  is  necessary  to  allow  for  the  heterogeneous  deviatoric  stress  field 
caused  in  fiber-reinforced  composites  by  the  hydrostatic  loading.  In  case 
of  anisotropic  materials  which  are  not  influenced  by  the  hydrostatic 
loading,  an  equation  similar  to  Equation  28  can  be  written  in  terms  of 
energies  of  distortion. 


Equation  28,  using  contracted  notation,  specialized  for  the  plane 
stress  condition  becomes 


m 


Ki[I0i<“i]  = 


I  ( i  =  1,2,6) 


(29) 


Using  the  results  of  tests  under  simple  load  conditions 


-m 


=  [/  °id‘i] 


(30) 


E  ill 


where 

and  shear  strains.  Combining  Equations  29  and  30,  the  following  equation 


"lu*  e2u  an(^  e6u  are  norma]  (tensile  or  compressive) 


is  obtained 


u 


CTj  d« 


yf  <rid*i] 

f  tu 


m 


=  I  Ci  =  1,2,6) 


(31) 


The  shape  of  the  failure  surface  (Equation  31)  in  the  strain  energy  space 
is  determined  by  the  shape  factor,  m;  for  m  =  2,  it  is  spherical  and  for 
m  =  1,  it  is  pyramidal. 
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SECTION  IV 
APPLICATIONS 

A  computer  program  embodying  the  concepts  of  Sections  II  and  III  was 
developed  to  determine  the  theoretical  responses  of  laminates  subjected 
to  in-plane  loads  to  compare  them  with  the  available  experimental  data. 

This  analytical-experimental  correlation  was  confined  to  data  for  the 
boron-epoxy  material  system.  The  basic  property  data  (results  of  simple 
tests)  of  the  boron-epoxy  material  system  needed  for  the  computer  program 
was  extracted  from  Reference  11.  In  the  application  of  the  technique, 
cubic  spline  interpolation  functions  (References  12  and  13)  were  used  to 
represent  the  basic  stress-strain  data  of  laminas.  This  functional  form 
of  stress-strain  curves  provides  a  much  better  representation  of  the 
stress-strain  curves  and  moduli  over  the  entire  range  of  the  curve  than 
the  two  or  three  parameter  mathematical  function  or  the  piecewise  linear 
approach  used  In  References  7  to  9. 

In  addition  to  a  suitable  representation  of  nonlinearity  of  the  basic 
stress-strain  curves,  the  criterion  defined  in  Section  III  was  used  to 
determine  the  ultimate  cumulative  response  of  laminates  subjected  to 
incremental  loads.  The  criterion  for  three  values  of  the  parameter  m, 
namely  m  =  1/2,  1,  and  2,  was  compared  with  some  of  the  other  failure 
theories  in  the  first  quadrant  of  the  stress  space  in  Figure  3.  The 
choice  of  the  value  of  the  shape  factor  m  for  a  material  system  depends 
upon  the  ability  of  the  criterion  to  match  the  experimental  results. 

At  present  there  is  no  reliable  experimental  data  of  boron-epoxy  unidirec¬ 
tional  laminas  under  biaxial  stress  states,  which  could  be  used  to  deter¬ 
mine  a  suitable  value  of  m.  However,  in  case  of  multidirectional  laminates 
a  better  corroboration  between  experimental  and  analytical  results  was 
observed  for  m  =  1.  For  this  reason  the  parameter  m  was  taken  to  be 
unity  for  the  study  reported  herein.  With  m  =  1  the  criterion  reduces 
to  the  maximum  strain  energy  theory  for  linearly  elastic  and  isotropic 
materials. 
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Figure  3.  Comparison  of  Various  Strength  Theories 
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The  hypothesized  criterion  assumes  strain  energies  under  longitudinal, 
transverse,  and  shear  loading  as  independent  parameters.  When  the  sum 
of  three  ratios  of  the  current  energy  level  to  the  available  one  equals 
unity  the  lamina  degrades  completely  (Equation  31).  The  three  terms  of 
Equation  31  at  the  time  of  failure  represent  contributions  to  degradation 
made  by  longitudinal,  transverse,  and  shear  stresses  acting  on  the  laminas. 
In  application  to  multidirectional  laminates,  the  criterion  is  satisfied 
plywise.  When  the  contribution  to  degradation  is  dominantly  transverse 
or  shear  and  characterized  by  matrix  failure  mode,  the  lamina  is  assumed 
to  be  capable  of  sustaining  longitudinal  loads. 

The  incremental  constitutive  law  in  conjunction  with  the  failure 
criterion  was  used  to  predict  the  response  of  various  boron-epoxy  lam¬ 
inates.  To  compare  the  predictions  with  the  available  experimental  data, 
we  examined  four  categories  of  symmetric  laminates: 

(a)  Unidirectional  laminates  Including  off-axis  coupons 

(b)  Angle  Ply  (^a)  laminates; 

(c)  (0/-a)  laminates;  and 

(d)  Multidirectional  laminates. 

1.  UNIDIRECTIONAL  LAMINATES 

At  present  the  technology  to  test  (0°)  laminates  under  biaxial 
stress  states  is  not  sufficiently  advanced  to  yield  reliable  data.  In 
the  absence  of  relevant  data,  no  analytical -experimental  corroboration 
for  these  laminates  was  feasible. 

2.  OFF-AXIS  COUPONS 

In  the  case  of  off-axis  coupons,  the  laminates  are  unidirectional 
but  the  load  axis  and  the  material  axes  do  not  coincide.  This  creates 
a  biaxial  stress  state  relative  to  the  material  axes.  Gripping  and  test 
techniques  being  employed  tend  to  Introduce  nonuniform  stress  states 
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within  the  test  section  of  the  coupon.  With  a  large  aspect  ratio 
(length/width  >  10),  however,  a  fairly  uniform  stress  field  develops  at 
the  midsection  of  the  coupon.  The  aspect  ratio  between  specimen  tabs 
used  in  the  experimental  program  (Reference  l'1)  from  which  this  data 
was  extracted  was  16.  The  specimens  were  0.5  inch  wide. 

In  developing  the  incremental  constitutive  relations  in  Section  II, 
we  assumed  that  the  lamina  was  orthotropic  in  its  behavior.  To  substan¬ 
tiate  this,  we  measured  the  stresses  and  strains  under  tensile  loadings 
of  off-axis  coupons  relative  to  the  material  axes;  these  values  are  plotted 
in  Figures  4  and  5,  along  with  simple  stress-strain  curves  obtained  from 
Reference  11.  The  resolved  stress-strain  curves  are  close  to  the  curves 
obtained  under  simple  load  conditions,  which  indicates  that  there  is 
practically  no  coupling  between  normal  and  shearing  deformations.  This 
leads  to  the  conclusion  that  the  lamina  does  tend  to  retain  its  orthotropic 
characteristics  to  failure.  In  the  absence  of  similar  data  under  com¬ 
pressive  loading,  we  assume  the  behavior  of  the  lamina  would  be  similar. 

Stress-strain  curves  of  15°,  30°,  45°,  and  60°,  and  75°  off-axis 
coupons  obtained  experimentally  (Reference  14)  and  using  analytical  tech¬ 
niques  are  shown  in  Figures  6  to  10.  Of  the  three  analytical  curves,  one 
is  based  upon  the  theory  presented  in  this  report  and  the  other  two  on 
the  concepts  of  References  7  and  8.  The  analytical  curve  based  on 
Reference  7  is  very  much  dependent  upon  the  size  of  the  increment  of  the 
load.  The  attempt  to  represent  curved  segments  of  the  basic  stress-strain 
curves  by  a  series  of  line  segments  introduces  errors  which  are  cumulative. 
Additionally  the  ultimate  cumulative  response,  according  to  this  technique, 
is  attained  when  any  of  the  strain  components  associated  with  the  material 
axes  reaches  the  limiting  value.  On  the  other  hand,  the  analytical  curve 
based  upon  Reference  8  is  not  dependent  upon  the  size  of  the  load  incre¬ 
ment,  but  it  makes  no  attempt  to  determine  the  ultimate  response.  This 
analysis  assumes  that  only  the  shear  stress-strain  behavior  of  the  lamina 
is  nonlinear.  Equation  1  can  approximate  the  curve  by  matching  the  slope 
at  the  origin  and  making  it  pass  through  some  other  convenient  point. 


18 


AFFDL-TR-73- 1 37 


Figure  4.  Lamina  Shear  Response 
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Figure  5.  Lamina  Transverse  Tension  Response 
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The  response  of  the  15°  off-axis  coupon  (Figure  6)  is  governed  mainly 
by  the  lamina  shear  stress-strain  response,  which  contributes  94.14%  to 
the  degradation  of  the  coupon  (Table  I  and  Figure  11).  The  stiffer  shear 
stress-strain  response  of  the  lamina  as  compared  to  that  of  the  15°  coupon 
(Figure  4)  is  reflected  in  the  difference  between  the  analytical  and  the 
experimental  behavior.  The  reasons  for  the  differences  in  the  analytical 
curves  were  discussed  earlier. 

With  increasing  orientation  of  fibers,  the  effect  of  the  transverse 
stress-strain  response  of  the  lamina  on  the  behavior  of  off-axis  coupons 
increases,  but  the  shear  effect  diminishes.  For  an  orientation  angle 
greater  than  60°,  the  transverse  effect  contributes  more  than  87.97%  to 
the  failure.  In  this  case  the  analytical  curves  based  upon  Reference  8, 
which  do  not  consider  the  nonlinearity  of  the  transverse  stress-strain 
response  of  the  lamina,  depart  more  from  the  experimental  behavior. 
Experimental  stress-strain  curves  agree  better  with  the  analytical  responses 
obtained  using  the  concepts  of  this  report  than  with  those  of  either 
Reference  7  or  Reference  8.  However,  experimental  and  predicted  strengths 
tend  to  diverge  with  increasing  orientation  especially  for  a>45°.  This 
may  well  be  due  to  such  factors  as  nonuniformity  of  stresses  in  the  test 
section,  the  width-effect  of  the  specimens  (0.5"  wide  specimens  were  used 
in  the  test  program),  residual  stresses^  etc . 

Theoretical  strengths  and  the  contributions  to  degradation  for 
various  orientations  of  fibers  in  off-axis  coupons  are  presented  in  Table 
I.  Percentages  contributions  to  degradation  of  off-axis  coupons  are 
shown  in  Figure  11.  The  theoretical  curves  in. Figure  11  indicate  that 

(a)  For  a  greater  than  7*,  matrix  failure  mode  dominates; 

(b)  For  a  =  10°  to  20°,  failure  is  essentially  shear  type; 

(c)  For  a  -  20  to  60°  failure  is  a  combination  of  tranverse 
and  shear  modes; 

(d)  For  a  greater  than  60°,  failure  is  a  transverse  one. 
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Figure  8.  Boron-Epoxy  45°  Off-Axis  Coupon 
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Figure  10.  Boron-Epoxy  75°  Off-Axis  Coupon 
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On  this  basis,  we  would  expect  that  for  a  =  10°  to  90°,  failure  planes 
would  be  parallel  to  the  fiber  oritentations  and  failures  would  corres¬ 
pond  to  matrix  failure  modes.  Figure  12  (References  14  and  15)  appears 
to  confirm  this  finding. 

3.  ANGLE  PLY  (+a°)  LAMINATES 


These  laminates  are  fabricated  from  unidirectional  laminas.  To 
e>at.e  their  analytical  behavior  to  the  basic  experimental  response  of 
laminas,  we  used  Equations  26  and  31.  Theoretical  failure  stresses  and 
contribution  to  degradation  of  0°,  ±15°,  ±30°,  ±45°,  and  (0790)  laminates 
subjected  to  various  stress  ratios  are  summarized  in  Tables  II  and  III. 

To  validate  these  results  would  require  an  experimental  corroboration; 
relevant  experimental  data  is  almost  nonexistant.  Reference  11  contains 
the  results  of  tests  conducted  on  (±30°),  (+45°),  and  (±60°)  laminates 
under  uniaxial  loads.  Some  experimental  data  from  biaxial  tests  on  (±45°) 
and  (0790°)  laminates  is  reported  in  Reference  16.  Analytical  stress- 
btrain  relations  of  (±45°),  (+30°),  and  (±60°)  laminates  under  uniaxial 
loads  were  obtained  using  the  present  theory  and  techniques  of  References 
7  and  9.  A  comparison  of  experimental  and  the  analytical  stress-strain 
responses  (Figures  13  to  15)  indicates  a  better  correlation  with  the 
present  theory.  Figure  16  indicates  that 


{a'  For  a  =  0°  to  30°,  longitudinal  contribution  to  degradation 
decreases,  shear  and  transverse  contributions  to  degradation  increase 
(transverse  stresses  being  compressive). 


J a  =  30°  1°  45°*  lonSitudinal  contribution  to  failure  con- 
bi?1  onCde?rlases?a "  COntr,bution  cont1™es  and  transverse 

(c)  For  a  =  45°,  shear  contribution  reaches  a  maximum  of  97.57% 
with  longitudinal  and  transverse  contributions  being  0.71%  and  1.72%. 


(d)  For  a  =  45°  to  60°,  shear  contribution  drops  and  transverse 
contribution  rises  sharply.  Transverse  contribution  is  due  to  tensile 
stresses. 

(e)  For  a  =  60°  to  90°,  the  failure  is  essentially  a  transverse 
failure  mode. 
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Figure  12.  Failed  Off-Axis  Coupons 
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Figure  13.  Response  of  (t45°)  Laminate  Subjected  to  Stress  a  Only 
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Figure  14.  Response  of  (+30°)  Laminate  Subjected  to  Stress  a  Only 
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Experimental  stress-strain  plots  for  (+45°)  laminates  subjected  to 
biaxial  loading  are  not  available,  so  no  comparison  is  possible.  A 
comparison  of  theoretical  and  experimental  strengths,  however,  for  (±45) 
laminates  under  biaxial  loads  is  shown  in  Figure  17.  The  results  do  not 
agree  too  well  with  experimental  data. 

Lacking  sufficient  and  reliable  experimental  data,  full  corroboration 
of  the  analytical  technique  of  Sections  II  and  III  as  applied  to  angle 
ply  laminates  is  not  feasible.  However,  some  insight  into  their  behavior 
may  be  obtained  from  the  analytical  results  summarized  in  Tables  II  and 
III.  Some  of  these  results  are  shown  in  Figure  18.  Figure  18  shows 
that  transverse  strength  tends  to  increase  with  increasing  angle  a  and 
longitudinal  strength  declines.  For  a  given  stress  ratio,  the  points  of 
maximum  strength  correspond  to  peaks  of  longitudinal  contribution  to 
degradation.  In  other  regions  either  transverse  or  shear  contribution 
to  degradation  dominates. 

4.  CROSS  PLY  (0790°)  LAMINATES 

A  special  case  of  angle  ply  laminates  is  (0790°)  laminate  with  equal 
number  of  0°  and  90°  plies.  Its  behavior  is  different  from  that  of  other 
laminates  discussed  so  far.  In  unidirectional  and  angle  ply  laminates, 
there  is  only  one  failure  state,  but  in  (0790°)  laminate,  such  a  state 
occurs  only  for  ax/ay  =1:1.  For  other  stress  ratios,  a Jo  ,  transverse 
degradation  of  0°  or  90°  plies  does  not  result  in  a  general  failure  of 
the  laminate.  Thus  (0°/90°)  laminate  has  two  failure  states,  the  initial 
failure  state  and  the  final  failure  state.  After  impairment  of  transverse 
load  carrying  capacity  of  0°  or  90°  plies,  the  affected  plies  tend  to 
unload  transversely.  The  manner  In  which  the  unloading  takes  place  is 
debatable.  In  the  analysis  reported  herein,  it  is  assumed  that  a  trans¬ 
verse  or  shear  Impairment  (l.e.,  on  the  occurrence  of  the  matrix  failure 
mode)  the  ply  or  plies  affected  by  such  impairment  unload  and  transfer 
transverse  and  shear  loads  to  the  unaffected  plies.  This  process  of 

load  transfer  is  assumed  to  continue  till  a  general  failure  of  the 
laminate  occurs. 
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Figure  17.  Strength  Envelope  for  (+45°)  Laminate 
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Figure  18.  Strength  Envelopes  for  (±a°)  Laminates 
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On  the  basis  of  the  above  assumption,  the  analytical  strength  envelope 
contributions  to  degradation  for  various  stress  ratios  and  stress-strain 
curves  were  obtained  (Figure  19  to  24)  and  compared  with  the  available 
experimental  data  (Reference  16).  In  addition  analytical  stress-strain 
curves  were  obtained  using  the  technique  of  Reference  7.  Though  the 
analytical-experimental  correlation  appears  to  be  reasonable,  there  are 
some  features  of  the  experimental  data  which  need  to  be  kept  in  mind  in 
evaluating  this  correlation.  These  features  are  the  differences  in  the 
properties  of  materials  used  in  the  experimental  and  analytical  programs, 
preload  effects,  and  the  reported  experimental  strain  data.  The  strains 
reported  are  the  readings  o*  the  three  elements  of  the  45°  strain  gage 

rosette  and  not  the  components  of  strain.  Hence  no  quantitative  compari¬ 
son  is  feasible. 

5.  (0°/+a°)s  LAMINATES 

In  this  category,  only  (0°/±60°)s  and  (0°/+45°)s  laminates  were 
studied.  Stresses  and  contributions  to  degradation  at  failure  for  these 
ami  nates  are  given  in  tables  IV  and  V.  Values  for  the  (0°/±60°) 
laminate  are  plotted  iri  Figures  25  and  26.  We  observed  that 

(a)  For  oy/ox<l ,  transverse  degradation  of  0°  plies  causes  initial 
failure,  but  the  laminate  is  still  capable  of  carrying  additional  load; 

(b)  For  ax/oy  *  1,  all  plies  fail  simultaneously,  resulting  in  the 
general  failure  of  the  laminate; 

(c)  For  ox/ay>l,  transverse  degradation  of  +60°  plies  causes  Initial 
failure,  which  is  followed  by  a  general  failure  when  0°  plies  degrade. 


No  biaxial  data  is  available  to  substantiate  the  analytical  results. 

A  comparison  of  experimental  results  obtained  from  testing  coupons  of 
(07+60°)s  laminate  under  uniaxial  loading  with  the  theoretical  results 
is  shown  in  Figures  27  and  28;  the  correlation  Is  very  good  for  the 

stress  ratio  ox/ay  ■  1:0,  but  it  is  far  from  acceptable  for  the  stress 
ratio  a  /a  «  0:1. 
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Figure  19.  Strength  Envelope  for  (0790°)  Laminate 
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Figure  22.  Stress-Strain  Response  of  (0790°)  Laminate 
Subjected  to  a  Stress  Ratio,  o  /o  “  1:2 
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Figure  27.  Stress-Strain  Response  of  (0°/+60°)  Laminate 
Subjected  to  Stress  only  s 
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Figures  29  and  30  show  the  analytical  strength  envelope  and  contri¬ 
bution  to  degradation  of  (0/±45°)s  laminate  subjected  to  various  stress 
ratios.  Since  no  experimental  data  is  available,  we  can  make  no  compari¬ 
son. 

6.  MULTIDIRECTIONAL  LAMINATES 

Innumerable  multidirectional  laminates  can  be  fabricated  from  uni¬ 
directional  plies.  For  this  study,  we  considered  only  (07+45790°) 
family  laminates.  There  are  three  laminates  in  this  category,  namely 
(07+45790° )$,  (0°2/+45°/90°)s,  and  (0°/±45°/0°/M°)s,  for  which  some 
experimental  data  is  available  for  comparison  with  the  analytical  results. 
Failure  stresses  and  contributions  to  degradation  for  the  above  three 
laminates  are  summarized  in  Tables  VI  to  VIII. 

Figures  31  to  36  show  the  strength  envelope,  contributions  to  degrada¬ 
tion,  and  stress-strain  curves  for  the  (0°/±45°/90°)  laminate.  The 
strength  envelope  is  symmetric  about  the  stress  ratio  of  1:1.  The  agree¬ 
ment  between  predicted  and  experimentally  determined  strengths  for  the 
stress  ratios  oy/ax  of  1:0  and  2:1  is  very  close  (Figure  31).  For  the 
stress  ratio  ax/ay  of  1:1  there  is  general  agreement  but  the  experimental 
values  have  considerable  scatter.  In  this  laminate,  the  initial  failure 
is  due  to  the  transverse  degradation  of  0°  plies  for  stress  ratios 
0<ax/ay<l.  The  final  failure  is  due  to  degradation  of  90°  plies  for 
stress  ratios  of  0<cx/o  <1/2  and  by  degradation  of  ±45°  plies  for  stress 
ratios  ax/oy  =  1/2  to  1.  For  the  stress  ratio  cx/oy  =  1:1,  all  plies 
fail  simultaneously;  the  qualitative  correlation  between  analytical  and 
experimental  stress-strain  curves  was  very  good,  as  shown  by  Figures  33  to 
36. 


In  the  case  of  the  (0°2/±45°/90°)s  laminate,  no  experimental  stress- 
strain  data  is  available.  For  this  reason  only  failure  stresses  (experi¬ 
mental  and  analytical)  are  compared  in  Figure  37.  The  agreement  between 
these  is  reasonable.  The  contributions  to  degradation  of  the  laminate 
for  various  stress  ratios  are  shown  in  Figure  38. 
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Figure  30.  Percentage  Contribution  to  Degradation  of  (0°/±45°) 
Laminate  Subjected  to  Various  Stress  Ratios 


:1gure  32.  Percentage  Contribution  to  Degradation 
Laminate  Subjected  to  Various  Stress  "L 
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Figure  33.  Stress-Strain  Response  of  (0°/t45°/90°)s  Laminate 
Subjected  to  Stress  o  Only 
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Figure  36.  Stress-Strain  Response  of  (0°/t45o/90o)  Laminate 
Subjected  to  a  Stress  Ratio,  0/0=  -1:2 
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Longitudinal;  T  =  Transverse;  S  =  Shear 


FAILURE  STRESSES  AND  PERCENTAGE  DEGRADATION  OF  (0°/±45°/90°)  LAMINATE 
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Figure  37.  Strength  Envelope  for  (0“/±45°/90°)s  Laminate 


Figure  38.  Percentage  Contribution  to  Degradation  of  ( 
Laminate  Subjected  to  Various  Stress  Ratios 
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Figures  39  and  40  show  the  theoretical  strength  envelope  and  contri¬ 
butions  to  degradation  of  (0o/±45o/0°/Mo)s  laminate.  The  experimental 
data  available  for  this  laminate  is  quite  meager  {Reference  16).  The 
comparison  of  this  with  analytical  data  {Table  IX)  indicates  good 

corroboration  for  stress  ratios  a  Jo  /t  of  1:0:0,  0:0:1,  1:2:1,  but  far 

x  y  xy 

from  satisfactory  correlation  for  the  stress  ratio  ax/ay  of  1:2. 

Strength  envelopes  of  three  (0°/±45°/90°)  family  laminates  are 
shown  in  Figure  41.  The  figure  indicates  that  the  laminate  increases  in 
strength  for  certain  stress  ratios  due  to  percentage  increase  of  0°  plies, 
but  decreases  in  strength  for  other  stress  ratios. 
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Figure  39.  Strength  Envelope  for  (07+4570°/M° )  Laminate 
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SECTION  V 

•DISCUSSION  OF  RESULTS  AND  CONCLUSIONS 

In  Section  II,  we  presented  a  technique  for  predicting  the  response 
of  composite  laminates.  The  incremental  constitutive  law  defining  the 
nonlinear  behavior  of  anisotropic  materials  is  specialized  for  ortho¬ 
tropic  laminas  under  plane  stress  condi tons.  Using  simplifying  assump¬ 
tions,  we  extended  the  use  of  the  incremental  constitutive  law  to  the 
analysis  of  multidirectional  laminates  subjected  to  inplane  loads.  In 
the  application  of  the  technique,  we  represented  the  nonlinearity  of 
basic  stress-strain  curves  of  laminas  by  the  use  of  cubic  spline  inter¬ 
polation  functions. 

The  incremental  constitutive  law  in  conjunction  with  the  failure 
criterion  was  used  to  predict  the  response  of  various  boron-epoxy 
laminates.  The  available  experimental  data  was  not  only  compared  with 
the  predictions  based  upon  the  present  theory,  but  also  with  theoretical 
results  obtained  using  techniques  of  Petit-Waddoups  and  Hahn.  In  the 
case  of  unidirectional  laminates  subjected  to  biaxial  stress  states,  we 
could  make  no  comparison  because  no  relevant  experimental  data  was 
available. 

The  study  of  other  laminates  revealed  that  the  boron -epoxy  material 
system  tends  to  retain  its  orthotropic  characteristics  to  failure.  Ur 
to  the  point  of  experimental  failure,  experimental  stress-strain  curves 
matched  analytical  curves  based  upon  present  theory  more  closely  than 
the  curves  using  the  techniques  of  References  7  to  9.  However,  in  the 
case  of  the  (0°/±60°)  laminate  subjected  to  the  stress  o,  considerable 

«  J 

deviation  between  experimental  and  analytical  results  was  observed.  The 
probable  cause  is  traceable  to  the  role  played  by  the  0°  plies  in  the 
laminate;  we  assumed  that  transverse  degradation  does  not  precipitate  a 
total  failure  since  the  plies  are  still  capable  of  sustaining  longitudi¬ 
nal  loads.  This  is  probably  not  realized,  and  the  0°  plies  may  lose 
their  total  effectiveness  on  transverse  impairment.  This  behavior  of 
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0°  plies  in  the  boron-epoxy  laminate  is  quite  different  from  that  observed 
in  a  glass-epoxy  laminate  (Reference  18). 

No  experimental  data  is  available  for  the  (07+45°)  laminate.  In 
light  of  the  behavior  of  (0°/+60°)s  laminate,  however,  we  expect  that 
the  transverse  degradation  of  0°  plies  would  precipitate  general  failure 
of  0°  plies  and  the  laminate. 


This  study  has  shown  that  the  incremental  constitutive  law  in  com¬ 
bination  with  the  failure  criterion  proposed  in  this  report  provides  a 
suitable  analytic  means  far  predicting  the  behavior  of  laminates  under 
general  states  of  stress  using  only  the  unidirectional  laminate  property 
data  obtained  under  simple  load  conditions.  Though  attention  was  con¬ 
fined  to  the  tension-tension  quadrant  of  the  stress  space,  the  technique 
is  general  and  holds  for  other  quadrants  as  we>l,  provided  no  buckling 
occurs.  It  has  also  demonstrated  that  a  need  exists  for  improving  the 
experimental  techniques  to  generate  good  data  to  clarify  fuzziness  that 
exists  between  analytic  and  experimental  responses. 
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